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The vacuum energy for a massless conformally coupled scalar field in a brane world set up, 
corresponding to de Sitter branes in a bulk anti-de Sitter spacetime, is calculated. We use the 
Euclidean version of the metric which can be conformally related to a metric similar in form to the 
^s^j ■ Einstein universe (S 4 x R). Employing ^-function regularisation we evaluate the one-loop effective 

potential and show that the vacuum energy is zero for the one brane and non-zero for the two brane 
configuration. We comment on the back-reaction of this Casimir energy and on the inclusion of a 
mass term or non-conformal coupling. 
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■ Many ideas and proposals within the circumference of what we call brane world cosmology (BWC) relate to the 
work of Randall and Sundrum (RS) JO who initially implemented their two brane model to try and solve the hierarchy 
problem. In a variant of this model (a single brane with positive tension 0) they also showed that the extra dimension 
could be of infinite extent. Generally speaking, the second model (and its extensions to include curved branes) has 
been of more interest to the cosmology community because gravity behaves essentially like 4-dimensional gravity 
localised near the brane . 

An interesting BWC scenario has been developed in ||, where it is possible to obtain inflation on a single positive 
tension brane solely due to the effect of a bulk gravitational scalar field. In this regard, the vacuum energy of the bulk 
scalar field could have some affect on the cosmological evolution of the brane, (depending on the size and sign of this 
quantity) and is an interesting subject in its own right. The importance of vacuum energy in brane world scenarios 
(BWS) has also been highlighted in || where it has been suggested that it may be possible to solve the hierarchy 
problem with two positive tension branes, if the back-reaction of the Casimir effect is included (also see ) . 

In this letter we evaluate the one-loop effective potential (leading to the Casimir energy) for one and two de Sitter 
. , brane configurations in a bulk 5-dimensional anti-de Sitter spacetime. A massless conformally coupled scalar field is 
investigated as this is the simplest case technically. In || a similar calculation has been considered for spherical and 
hyperbolic one brane configurations. However, we consider a spherical two brane configuration and for the one brane 
case the results agree. Quite extensive work has been done for flat branes and we refer the reader to the literature, 
see P-[l4[ and the references therein. 

Calculations of the Casimir energy are non-intuitive in how the overall sign depends on the topology and/or 
dimension of the spacetime |l5| . One has to perform the calculation. To obtain the result, we employ the ^-function 
regularisation method Jl6| |. The beauty of (^-function regularisation (or some would say its weakness) is that it 
automatically removes the divergences of the theory leaving a finite result. An interesting discussion about different 
analytic continuations in £-function regularisation can be found in Q. 

To perform our calculation we use the Euclideanised form of the metric suitable for de Sitter branes . It can be 
written in conformal form as 

ds 2 = a 2 (z)(dz 2 + d X 2 + sin 2 X dfi? 3) ); a(z) = -— -r-rr , (1) 

' smh(z + \z\) 

where t = (— 6/A5) 1 / 2 is the anti-de Sitter radius, dfl 2 ^ is the metric on the 3-sphere and the positive tension brane 
is placed at z = 0. Here, we place the other negative tension brane at \z\ = L (the one brane configuration is given 
by L — > 00). The brane tensions in this case are |lg| , |l9[ | 
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and the non-dimensional length L is given in terms of the physical length r as 

dr 
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r + dr 
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The metric is conformal to the Euclidean version of an Einstein universe (5 4 x R) with unit radius. The topology 
of this metric corresponds to two concentric hyperspheres of equal radius on a cylinder of length L (the conformal 
scaling has shrunk the radii of the two hyperspheres to equal radius). 

Because we are dealing with a conformally coupled scalar field we must 

find the eigenvalues of the Klein-Gordon equation given by (e.g., see M) 



(-9 z 2 -A( 4 )+e 5 i? (4) )0 = A 2 0, 
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A (4) 



cj) = X 



(4) 



where £5 = 3/16, R 1 - 4 ^ = 12 and imposing Z2 symmetry on the branes implies Neumann boundary conditions giving 
the discrete sum between the brane boundaries, for a conformally coupled scalar field p0| . In general, as pointed out 
in real the eigenfunctions will satisfy mixed boundary conditions and the spectrum of —dp' will not be easily obtained. 



For the operator — A^ 4 ) +9/4, the eigenvalues are (for Gegenbauer polynomials see pi].) 
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X 2 m = (m + 3/2) 2 , 



d(m) 



■(m + 2)(m + 3/2)(m+l), 



where dim) is the degeneracy. Thus the eigenvalues for the Klein-Gordon operator are 



(^) 2 + (- + 3/2) 2 



Employing the Euclidean effective action, we define the function 

00 

m,n— 

and then the one-loop effective action, W\, is related to g(s) by 



Wi 



(5) 



(6) 



(7) 



(8) 



Note, there is no renormalisation scale, fj,, in the function g(s) because the eigenvalues are dimensionless (due to the 
conformal rescaling). 

We first begin with the simpler case of the vacuum energy for the Einstein universe S x R. This is equivalent to 
finding the zero point energy of the equation 



a?V - a (4) ^ + — i? (4 V = 0. 
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The eigenvalues and degeneracy are still given by Eq. (0). It is then straightforward to show that 
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E = - g(-l/2) = - { C(2s - 3, 3/2) - 7 C(2.s - 1, 3/2) 



= 0, 



(9) 



(10) 



-1/2 



where £(a, b) is the generalised (or Hurwitz) ^-function and is related to the Bernoulli polynomial B n {x) by plj 

C(-n,a) = — — . (11) 

n + 1 



It is easy to verify that E is zero, using Eq. (11), as should be expected because in an odd number of dimensions 
there is no conformal anomaly EQ] . This calculation is similar to the one by Ford [E2[ for S 3 x R, where an exponential 



cut off was used in the mode sum rather than C-function regularisation (for S 2 x R see 23 1). 

We now evaluate the effective potential for a one brane configuration, where L — > 00 in Eq. (^) and the discrete 
sum over n is now an integral. The function g(s) becomes 
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9i(s) - — / dk V d(m) (fc 2 + (m + 3/2) 2 p , (12) 

where the factor of 2 is because there are two copies of the bulk space on either side of the brane. For large s we can 
interchange the order of the sum and the integral and perform the k integration, leading to 

= ^ £ ^P (C(2- - 4, 3/2) - \((2s - 2, 3/2)) , (13) 

where in the second step we have used simple algebra to rewrite the equation in terms of generalised ^-functions. For 
s = it is clear that g\ (0) = because 

' .s + 7 s 2 + 0(s 3 ), (14) 



where 7 is Euler's constant. Thus, 

5i(0) = ^ r (-V2) (c(-4,3/2) - Jc(-2,3/2)) = 0, (15) 

which is zero because the combination of £- functions in Eq. ( p"5| ) cancel, as can be verified by employing Eq. (pj]). 
Therefore, the one-loop effective potential (and hence the vacuum energy) for a conformally coupled scalar field is 
zero for the one brane case and does not affect BWC models, such as in ||. This is consistent with the result found in 
pL Although, one might naively expect imposing a boundary on the brane similar to that for a half Einstein universe 
will give a non-zero contribution. 

Now we come to the main task at hand. The two brane configuration has two infinite summations, which can be 
expressed in terms of generalised functions. In fact, there are many subtle issues regarding the correct analytic 
continuation of such functions upon interchange of the order of the summations (see |lq ] for a detailed discussion) . 
The function g(s) becomes 



00 00 



.9200= E E d ( m K cV + ( ra + 3 / 2 ) 2 

n— — oo m— 

1 00 00 poo 

= Y(s) E E d ( TO )/ dtt s - l C M-t[c 2 n 2 + {m + 3/2f]}, (16) 

* n=-oo m=0 



where c = tt/ L and in the second step we have made use of the Mellin transform. The sum over n is from —00 to 
00 because imposing orbifold boundary conditions essentially doubles the modes. This formula is almost identical 
in form to the 2-dimensional Epstein-Hurwitz £- function studied by Elizalde in |24|] (see also |l(|), apart from the 
degeneracy factor dim). This work has many interesting results and the one we employ is 



00 1 — I — 00 

E e - c2 " 2 = V? + 2 V?^ exph7rV/c2] - 



(17) 



92(s) = +W72- tv„/ ( «2s - 4, 3/2) - -C(2s - 2, 3/2) 



This is non other than the unimodular transformation formula of Jacobi's 63 function [ pl| (the Poisson resummation 
formula). Substitution of the above equation into ([jii]) allows us to interchange the order of the summations. After 
integrating with respect to t and performing some simple algebra we obtain 

1 IT T(s- 1/2) ( \ fn A 1 
3V ^ W) 

+ ^rrc- s - 1/2 V d(m) n s - 1 / 2 (m + 3/2)- s+1 / 2 ^ s _ 1/2 [2^c- 1 n(m + 3/2)]. (18) 
T(s) f-^ 

n— 1, m=0 

It is easy to verify that 32(0) = 0, using Eq. (|l4|). Furthermore, taking the derivative of 172 (s) with respect to s (and 
leaving only terms that remain independent of T(s) near s = 0, see Eq. dl4)) gives 



92(0) = +^^(-1/2) ^(-4, 3/2) - jC(-2, 3/2)) 

OO 

+ 4c_1/2 E d(m)n- 1 / 2 (m + 3/2) 1 / 2 K_ 1/2 [2irc- 1 n(m + 3/2)}, (19) 



n—1, 771—0 
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where a prime denotes differentiation with respect to s. Interestingly, the first term is exactly the contribution from 
the one brane configuration Eq. (|l^), which is zero. This is similar to finite temperature field theory where the n = 
mode gives the zero temperature contribution (this can be clearly seen by using heat kernel methods psj] , see also 
H). Of course, this kind of behaviour was expected due to the topology of the conformally re-scaled metric, Eq. (Q). 
The second term is non-zero and depends on the value of c = ir/L. Note that K_i/ 2 (z) — ^JirJ(2z) e~ z . Therefore, 
the one-loop effective potential per unit volume is given by, after some simplification, 

v Wi A(L) 

1 2LVol(S 4 ) 8tt 2 L' 1 ' 

where W\ is the effective action for the conformally re-scaled field and 

(m + 2)(m + 3/2)(m+ 1) 



A(L) = 2^ '—^ exp[-2Ln(m + 3/2)] 

n—1, m—0 
oo 

^2 (m + 2)(m + 3/2) (m + 1) ln(l — exp[— L(2m + 3)]) . (21) 



m=0 



The function A(L) is generally a complicated function of L, but its qualitative behaviour can be deduced analytically. 
First, it is easy to see that A{L) > and dA(L)/dL < for L > 0. Note, that we have normalized the radius of S 4 to 
unity. Therefore, the limit L — > corresponds to infinitely large 4-spheres on the boundaries, i.e., flat space. In the 
limit L <C 1, we may approximate A(L) as 



/>oo 

A(L) ~ - / dmm 3 ln(l - exp[-2mi]) 
Jo 



which gives 



Using the relation, 



4(2£) 4 J e x -l 4(2L) 4 



c- ( -2„) =,-!)■ rl2 ";g-;i 2 r" ^.,, H , ^ 



we find Eq. ( p3| ) exactly agrees with the flat space Casimir energy for a conformal scalar JTo[ ] . It may be worth 
mentioning that V\ given by Eq. ( p3| ) is also equal to the free energy of a massless boson at temperature T = 1/(2L), 
except for the factor of 2 that comes from the mode doubling. In the limit L ^> 1, A(L) behaves as 

oo 

A(L) «^(m + 2)(m + 3/2) (m + 1) exp[-£(2m + 3)] - 3exp[-3i] . (25) 

Thus, the Casimir energy for de Sitter boundaries decreases much more rapidly than that for flat space boundaries 
as the non-dimensional length L increases. We have checked our results for the sum A{L) on different numerical 
packages and find that, for example, when L = 0.5, 1, and 2, A(L) = 5.90, 0.31 and 0.008, respectively. 

From the symmetries of the metric, Eq. (Q), and because the stress energy tensor for a conformally coupled scalar 
field is traceless in odd dimensions (no conformal anomaly, see [ppf) the Casimir energy density, p, is equal to the 
effective potential, V\. Also, using the fact that the conformally re-scaled stress energy tensor is given in terms of the 
original one by (T^ v ) g = a-° (T^„), see we have 

(T» v ) g = r 5 sinh 5 (z + \z\) (T» v ) . (26) 
This then implies that the physical energy density p g is given by 

sinh 5 (z + \z\) 

Ps = Stt^l A{L) ■ (2?) 

We now discuss the implications of these analyses on the issues raised in the introduction. 

The back-reaction of the Casimir energy plays an important role in BWS PJ^l, including inflating branes. However, 
in p] the back-reaction for the de Sitter brane model only considered the vacuum energy from Minkowski branes. For 
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flat branes this only depends on the difference between the number of bosons and fermions (Nr), where Ni, = Nt 
gives the usual RS model and Nt > Nf gives two positive tension branes, see ||. Thus, we should also consider other 
fields such as the fermion. However, in contrast to Minkowski branes, the eigenvalues in the mode sum for two de 
Sitter branes are different for bosons and fermions, see p2|| for S 3 x R, implying the vacuum energy will not be of 
equal magnitude, but opposite in sign. Therefore, we can infer, from the previous sentences, that the bulk will not 
remain pure anti-de Sitter under the back-reaction (for de Sitter branes) even for Nt — Nf. 

In the context of BWC models, such as in we must consider the more general case of a massive or non-conformal 
scalar field. Adding a mass term will break conformal invariance. Altering the conformal coupling constant £5 not 
only breaks conformal invariance, but also modifies the boundary conditions. As we mentioned earlier (also see ||), 
the boundary conditions will become mixed. One can of course consider a small perturbation of £5 or a small mass 
term. Such cases may be treated analytically. Otherwise, the eigenspectrum of —&£ needs to be solved numerically. 
Nevertheless, for large n the spectrum will return to irn/L and we expect to need only evaluate this correction to the 
Casimir energy up to, say, n « 10. 

In forthcoming work p5| we will present a more detailed discussion of the ^-function regularisation procedure. 
There, its relation to the heat kernel will be made explicit. In connection with the generalisation to non-conformal 
cases, we mention the work done by Dowker and Apps p6| for generalised cylinders (up to 4-dimensions), and the 
work by Garriga, Pujolas and Tanaka [|l0|j27| l for a 5-dimensional spacetime bounded by flat branes. They stress 
the importance of the cocycle function. Also, we will present the results for a spinor field and an analysis of the 
introduction of a small mass term (or non-conformal coupling) into the wave equation. Other avenues of research for 
de Sitter branes include the general non-conformal case, finite temperature effects (see |28| for the RS model and also 
the comments made in p9|), boundaries on the brane, and the possibility of stabilising the radion field. 
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